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A colored weak singlet scalar state with hypercharge 4/3 is one of the possible candidates for the 

explanation of the unexpectedly large forward-backward asymmetry in ti production as measured 

by the CDF and D0 experiments. We investigate the role of this state in a plethora of flavor 

changing neutral current processes and precision observables of down-quarks and charged leptons. 

Our analysis includes tree- and loop-level mediated observables in the K and B systems, the charged 

'"H lepton sector, as well as the Z ^ bb decay width. We perform a global fit of the relevant scalar 

^~~^ couplings. This approach can explain the {g — 2)^ anomaly while tensions among the CP violating 

^^ observables in the quark sector, most notably the nonstandard CP phase (and width difference) in 

Cn the Bs system cannot be fully relaxed. The results are interpreted in a class of grand unified models 

^ which allow for a light colored scalar with a mass below 1 TeV. We find that the renormalizable 

O 5*17(5) scenario is not compatible with our global fit, while in the SO(IO) case the viability requires 

^^ the presence of both the 126- and 120-dimensional representations. 

0\ PACS numbers: 14.80.Sv, 12.15.Ji, 12.10.-g, 12.10.Kt 
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I. INTRODUCTION 



Recent CDF and D0 results on the forward-backward asymmetry (FBA) in top quark pair production have 

attracted a lot of attention and a number of proposals have been made in order to explain all the relevant observables 

(for a recent review see [1]). Among these, a colored weak singlet scalar with charge 4/3 (3, 1, 4/3), if exchanged in the 

vj M-channel, can well accommodate most of the present measurements [2, 3] (see however also [4, 5]). Motivated by the 

J^ success of this proposal [2] we have systematically investigated the role of such state in charm and top quark physics [6] . 

Qv Constraints on the relevant couplings to up-type quarks come from observables related to D" — D'^ oscillations, as well 

^y^ as di-jet and single top production measurements at the Tevatron and the LHC. In turn, we were able to predict the 

l/~) expected rates of flavor changing neutral current (FCNC) mediated top quark and charmed meson decays, generated 

1^ by the presence of the new colored scalar. Due to the possibility to accommodate such states within SU{5) grand 

(^ unified theories (GUTs) that contain 5- and 45-diniensional Higgs representations, we have determined the resulting 

T— I textures of the up-quark mass matrix at the GUT scale. The particular SU{5) model we advocated has an appealing 

y—i feature of correlating the presence of light colored scalars stemming from the 45-dimensional Higgs representation with 

IJ bounds on the partial proton lifetimes. Namely, the aforementioned representation contains among other states two 

• ^H colored scalars — (3, 1, 4/3) and (8, 2, 1/2) — whose masses are below or of the order of 1 TeV when partial proton decay 

^^ lifetimes are predicted to be at or slightly above the current experimental bounds. The most common renormalizable 

H models based on 5*0(10) framework [7, 8], on the other hand, usually rely on inclusion of 120- and 126-dimensional 

scalar representations to generate fermion masses. As it turns out, both of these contain a colored weak singlet 

(3, 1,4/3) state that, if light, could accommodate the ti production observables [2, 9]. Our analysis of its couplings 

could thus help in establishing the viability of such 50(10) scenarios as we demonstrate later. 

During the last decades, rare processes involving down- type quarks and charged leptons have played an important 
role in revealing possible signs of new physics (NP) at low energies. A prominent example is the anomalous magnetic 
moment of the muon, whose most precise experimental measurement [10] deviates from theoretical predictions within 
the SM [11] by about three standard deviations. Also most recently, the CDF and D0 experiments [12-15] have 
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reported indications of a large CP-violating phase entering the Bg — Bg mixing amphtude, which cannot be accom- 
modated within the standard Cabibbo-Kobayashi-Maskawa (CKM) framework.^ Furthermore, the experimental data 
seem to prefer the decay widths' difference between the Bg eigenstates AF^ larger than predicted within the SM. 
Any NP addressing this discrepancy would have to contribute to the absorptive part of the mixing amplitude [20]. 
Finally, the recent measurements of the leptonic B ^ tv branching ratio induce a 2.9 ct tension on the global fit to 
the CKM unitarity triangle which may be ameliorated via a small NP contribution to the Bd — Bd mixing phase [19]. 

Motivated by the interesting role the scalar (3, 1, 4/3) state might play in down-quark and charged lepton physics, we 
systematically investigate its contributions to the down-quark and charged lepton physical observables. We consider 
observables affected already at the tree-level, as for example KL,s{Bs,d) -^ (-^f- ~ , lepton flavor violating (LFV) tau 
decays and jjl — e conversion in nuclei, and also loop suppressed contributions to e/f , Am^f , Bd^s — Bd.s oscillation 
observables, lepton anomalous magnetic moments, LFV radiative tau and muon decays, and the Z ^ bb decay width. 
These constraints can be consistently implemented within the framework of SU{5) and 5*0(10) GUTs which rely on 
tree-level generation of charged fermion masses. 

The outline of this article is as follows. In Sec. II we define the leptoquark couplings of the scalar (3, 1, 4/3) to SM 
fermions. In Sec. Ill we study the effects of this state on low energy precision observables and perform a global fit of 
its couplings in Sec. IV. Resulting values of couplings are then reanalyzed in Sec. V, where we study how they relate 
to the mechanism of mass generation in GUT scenarios and derive bounds on the vacuum expectation value of the 
45-dimensional Higgs representation in the SU{5) case. Finally we conclude in Sec. VI. 

II. ELECTROWEAK SCALE FRAMEWORK 

We consider a color triplet, weak singlet scalar with charge 4/3 

A -(3, 1,4/3), (1) 

which can couple to the right-handed fermions of the SM via the interaction Lagrangian 

Ca = Y,,hPLd%^''* + ^-^ea^>cU^aPLu'^b^' + h.C. , (2) 

where Pl.r — (IT75) and the totally antisymmetric tensor €abc is defined with 6123 = 1. The above interaction terms 
exhaust all possibilities of renormalizable couplings between SM fermions and the A scalar. As mentioned in Sec. I, 
the diquark couplings g.^- of A to up-quarks of different generations UiUj can play an important role in top and charm 
physics. The leptoquark nature of A, on the other hand, is parameterized by couplings Yij to charged leptons and 
down-quarks, (.idj. If, and only if both g and Y are present, baryon [B) and lepton (L) numbers are violated while their 
combination B — L is conserved. Proton lifetime is protected in this general case by antisymmetric color contraction 
between A and two up-quarks. It implies fiavor-space antisymmetric coupling gij{= —9ji) which prohibits the proton 
from decaying via dimension-6 effective operators mediated by A. To comply with the ti production parameters, the 
mass of A should be below 1 TeV, preferably around 400 GeV. This setup is natural in a theoretically well-motivated 
class of grand unified models. A realistic GUT context of (2) and the resulting model building constraints will be 
presented in Sec. V. 

In general y is a complex matrix acting on charged-lepton and down-quark flavor indices £RY{dji)^ . From the 
right-hand side of Y one may redefine the quark fields using the global B symmetry transformation, which since it is 
broken, has a side effect of globally rephasing diquark couplings g while leaving the mass and CKM matrices invariant. 
This is not a worry at this point since current experimental constraints from top quark physics and D^ — D^ mixing 
observables cannot probe the overall phase of 5 [6]. One can redefine lepton fields in an analogous manner. However, 
of the two independent phases used to redefine quark and lepton fields only their sum is physical, while their difference 
corresponds to B ~ L, a conserved quantum number. As a result, freedom remains to choose one phase in Y. Strictly 
speaking, from the phenomenological point of view in Sec. Ill where we do not consider observables sensitive to lepton 
mixing, we could have rephased charged lepton flavors independently. This would allow us to rephase each row of Y 
separately. We restrain however from using this freedom which would result in the leptonic phase convention being 
"gauged" according to Y instead of to the standard form of the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix. 



^ The recent D0 measurement of an anomalous like-sign dimuon charge asymmetry [16, 17] is consistent with the hypothesis of NP 
contributing only in B^ g — B^ ^ mixing [18, 19]. 



III. LEPTOQUARK PROBING OBSERVABLES 

The leptoquark couplings endow the scalar A with a potential to cause large effects in (flavor changing) neutral 
current processes of down quarks and charged leptons (see [21] for a recent analysis of scalar leptoquark constraints 
from K and B sectors). The couplings Yij of Eq. (2) must therefore pass constraints from plethora of precisely 
measured or bounded low energy observables. In this section we make predictions of the observables most sensitive to 
effects stemming from A and compare them to current experimental values. For each observable we state an effective 
error of the constraint, which is, as will be explained in the following, a total combined theoretical and experimental 
uncertainty. In order to confront this model with experimental data in a quantitative manner we wrap up the analysis 
with a global fit of all the 9 entries of Y in Sec. IV. 

A. SM theoretical inputs 

Most observed flavor phenomena are well described within the SM and thus the allowed size of NP contributions 
crucially depends on reliable estimates of SM parameters. In the presence of NP virtual contributions to quark flavor 
observables the extraction of the CKM matrix becomes more involved, since some observables used in the conventional 
CKM fits receive contributions from both the SM and NP amplitudes. As we want to treat the SM contributions as a 
theoretical background, it is imperative to calibrate the CKM matrix exclusively on SM tree-level observables, which 
are largely insensitive to virtual A contributions. 

Thus we employ the results of a simple CKM fit to tree-level observables.^ These are the measurements of the first 
and the second row CKM element moduli from super-allowed /3 decays, leptonic and semileptonic meson decays, as 
well as the extraction of the CP phase angle 7 from tree-dominated B decays [23] 

/0.97425(22) 0.2252(9) 3.89(44) x lO^^X 
JVckmI = 0.23(11) 1.023(36) 4.06(13) X 10-2 , 7 = 73(+^2)°. (3) 

In particular, the value of \Vub\ is an average of exclusive and inclusive semileptonic B decay analyses. We explore 
the impact of the branching ratio oi B ^ tv on the CKM fit in Sec. IV B. Note that we cannot use direct |Vt;,| 
determination from single top production measurements, since these may be affected by A contributions [6]. By 
fitting constraints (3) to the Wolfenstein expanded CKM matrix up to order A^, we find values in agreement with [22] 

A = 0.22538(65) , (4) 

A = 0.799(26) , 
p = 0.124(70) , 
7? = 0.407(52) , 



while we also extract the correlation matrix between the fit parameters 

(5) 



1 
-0.178 1 

-0.00517 -0.0553 1 

-0.0226 -0.242 -0.198 l/ 



In addition, we use the top quark pole mass of mj = 173.3 GeV [24], and the MS bottom and charm quark masses 
mbi'iTT'b) = 4.2 GeV, and mcimc) = 1.29 GeV [23]. Observable-specific numerical inputs will be stated where needed. 

B. Tree-level constraints 

We first focus our attention on observables which receive possible A contributions already at the tree-level and 
thus represent potentially most severe constraints on the Y matrix. The relevant effective Lagrangian for processes 



^ We do not use available results in the literature since they do not provide correlations among the parameters (e.g. UTFit tree-level 
fit [22]). 



involving charged lepton and down-quark pairs results from integrating out the A at tree level. After applying Fierz 
identities we recover the LFV and quark FCNC interaction terms among the right-handed leptons and quarks 

^%S,^,,,t = -^i^^^'P^^'')(dn,Pndd . (6) 

The corresponding leptonic (LFV) decay width of a neutral pseudoscalar meson P{didj) -^ ^a^t ^^ given by 

rp(...-,)^,,.+ = ^f^ "^ [^l + <- iK - rnlf] [(1 - (m. 4- fn,f){l - (m. - m,f)] '" , (7) 

where mp is the decaying meson mass, its decay constant is defined as customary for light neutral mesons {tt'^,K^), 
(O I dj^'^^^di I P{p)) — ip^fp/V^, while the hatted masses of leptons are rha^b = 'ma,b/'mp- We study the particularly 
interesting decay modes below. 

1. Kl — >■ lJ-~ fJ''^ , e+e~, fi'^e^ 

While the decay Kl — > 1^1^^ has been measured with great precision {B = (6.84±0.11) x 10~^ [23]), the presence of 
long-distance intermediate states Kl -^ 7*7* -^ A*^/^^ precludes similarly reliable SM predictions for this observable. 
We use a conservative estimate for the pure short distance branching fraction B'^ < 2.5 x 10~^, obtained using 
dispersive techniques [25, 26], as a Icr upper bound. Since the SM short distance contribution SgM(sD) ~ 0-9 ^ 10~^ 
is much smaller, we can neglect it and keep only the A-mediated amplitude. For the decay width K]^ — )■ n^n~, CP 
violation in K ~ K mixing is irrelevant and we treat Kl as a pure CP-odd state. Contributions of both K'^ and K'^ 
amplitudes are to be taken into account using Eq. (7) by replacing YajY^i* — >■ V2Re(Fps^pd*)- The decay width, 
mediated by CP conserving combination of couplings Y, is then 
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K.^.-,^ - — ^ — ^™^Vl-4-^ (8) 



A 

Using the lattice value of the kaon decay constant Jk = 156.0 MeV [27], the numerical result for the la upper bound 
is 

[Re(V„va]=<2,7xl„-»(j,^)\ ,9) 

In the di-electron mode K^ -^ e+e~ the experimental measurement of S = (Sj^^) x 10~^^ [23] agrees well with the 
long-distance dominated SM estimate oiB^ = (9 ±0.5) x 10~^^ [28]. The A contribution to this decay mode can be 
obtained from (8) by replacing ^ with e everywhere. Saturating the experimental uncertainty leads to the following 
1 a constraint 

[Re(y.r;,)]^<2.5xio-(^^^)\ (lo) 

A much stronger upper bound of the LFV decays B{Kl -^ M^e^) < 4.7 x 10^^^ at 90% confidence level (C.L.) has 
been set in [23]. The corresponding form of Eq. (7) is obtained by adding first coherently the flavor components of 
Kl and then summing over the widths of the two oppositely-charged final states. The result, with rUe set to zero, is 

r _t _ lYf^sY*^ + Yf,dY*J ruj^fj^ 2 r-, ^2]2 ,,.x 



5127r 



'A 



and implies a 1 ct bound 



\Y,sY:, + Y,,Y:f < 1.2 X 10-" (40^)' ■ (12) 



decay mode 90% C.L. exp. bound on B la upper bound in units (7nA/400 GeV)* 
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2.8 X 10"^ 
2.0 X lO"'' 



\Y,^Y:f < 4.4 
\Y,,bY;a\^ < 5.0 X 10- 
\Y^tY:f < 1.3 X 10- 

\Y,iY:f < 10.1 
KbY^f < 1.1 X 10- 



\YebY;d\ +\Y^tY:d\^ < 1.6x10- 
\Y^bY:d\^ + \Y^bY;/ < 2.2x10- 

inf,i;di' + \Y,bY:f < 2.7 x 10- 



J'^eb-'^lis ~r \if_ib^e 



< 3.4 X 10- 



Table I: Limits on Y couplings coming from upper bounds of lepton flavor conserving and violating B, 
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ri+ decays [23, 31]. 



2. Ks 
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Since Ks is approximately CP-even and is decaying to a CP-odd final state this decay mode is sensitive to the 
imaginary parts of y. In the muonic channel, the best limit still conies from the early seventies with B{Ks — >■ fj.'^ fi~) < 
3.2 X 10"'^ at 90 % C.L. [23], while the best upper bound on the branching fraction B{Ks -> e+e") < 9 x 10"^ at 90 % 
C.L. was more recently set by the KLOE experiment [29]. Both are still far above the SM expectations, whose long 
distance effects through Ks -^ 7*7* -^ e'e+ifi+^i-) reach 8 x 10"'3(2 x 10"^) x BiKs -^ 77) ~ 10""(10"") [23, 30]. 
These observables thus present clean probes of CP violating effects in the effective Lagrangian (6), through the decay 
widths 
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(13) 



(14) 



The resulting bounds, although diluted by helicity suppression in the electron mode and the short lifetime of Ks, are 
important constraints to be fulfilled by the following combinations of couplings at Icr C.L. 



[im(i;,y;,)]2 < 0.13 



mA 



400 GeV 



[Im(i;,r;,)]' < 1.1 X 10-4 ( 



TOA 



-V 

400 GeV/ 



(15) 
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In the SM these FCNC processes suffer additional helicity-suppression {m'j/m'^) leading to branching fractions of 
the modes with electrons, which are negligibly small compared to the current sensitivities of experiments, as given by 
the 90% C.L. upper bounds on B{Bd -^ e"e+) < 8.3 x 10"^ and B{Bs -^ 6^6+) < 2.8 x 10"^ [23]. In the dimuon 
channel the SM predictions for the branching fractions — of order ^ 10"^*^ (10^^) for Bd {Bg) decays — are closer to 
but still an order of magnitude below current experimental 90 % C.L. upper bounds 4.2 x 10"^ (1.2 x 10"®) [31]. Even 
in the case of Bd -^ t^t^ where the helicity suppression is the least severe, the SM prediction oi B ^ 10"^ [32] is 
far below the current experimental reach of 4.1 x 10"^ at 90% C.L. [33]. Consequently we do not need to consider 
pure SM or interference terms between SM and A-mediated amplitudes and focus our attention only to the pure A 
contributions. 

In Eq. (7) we substitute fp — > y/2fB^,^-. in order to conform with the standard normalization of heavy pseudoscalar 
decay constants. The lepton flavor conserving decay widths then read 
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while the rates of LFV decays are, e.g., for the ^t final state 
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(18) 



For the other dilepton LFV decays one should adapt the lepton indices of Y and replace rhr with the mass of the 
heaviest lepton in the final state. For the decay constants we use the central values of recent lattice QCD averages [27]: 
fsd = 193 MeV, fs^ = 239 MeV. The compilation of experimental upper bounds and their resulting interpretation as 
constraints on Y are given in Tab. I. 

4. B ^ Xs£+r 

EfFective Lagrangian (6) also contributes to the non-helicity suppressed b — >■ si~^£~ transitions. In particular it 
contributes to the Cg/ and Cio' Wilson coefficients of the effective weak Hamiltonian as defined in [34] . Following this 
reference, we write 
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where A^ -^ = VtbV^g, and A only contributes to 
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at the tree-level with the weak-scale Wilson coefficients 
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^■KYgYl, 
^GF\^t^OLm\ 



(21) 



The running of the Og, j^q, operators from the weak matching scale to the 6-quark mass scale is dominated by 
electroweak effects [35] and can be safely neglected for our purpose. At present the most sensitive observable is 
the inclusive decay width of B — > Xs£+£~, where I = e,/i, integrated in the dilepton invariant mass range of 
mf+£- = y^{pe+ +Pi-Y € [1,6] GeV. The corresponding branching fraction (S(i_6)ggv) is known in the SM to 10% 
accuracy and can be written in presence of Cg/^io' contributions as [34] 
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(15.86 ± 1.51) -0.049Re(C^,) + 0.061Re(Cfo,)+0-534|C^,|2-f0.543|Cfo,|2] x 10-^ (22) 



The experimental measurements of this quantity by the BaBar 

?cxp 
■^(l-6)GoV 



muon and electron flavors, yield 'B^^g^Q^y — (1.60 ± 0.5) x 10 ^ [38]. 



and Belle [37] experiments, averaged over the 



5. B ^ Tv£+£'- and B -^ Kl+i'' 

The exclusive B — >■ 7r£+^^ mode, where £ = iJ,,e, is severely CKM suppressed in the SM leading to branching ratio 
predictions which are well below the present experimental bound B{B^ — > ir'^t't^) < 4.9 x 10~^ @90% C.L. [23]. In 
addition, several LFV B^ -^ 7r+(A'+)£+£'^ modes have also been searched for at the S-factories and we compile the 
present bounds in Tab. II. 

The computation of the A contributions to these exclusive rare semileptonic B decays requires the knowledge of 
the relevant hadronic (7r|j7^|i?) and {K\J'^\B) matrix elements, where J^^ — bjfj^Pnq is the relevant quark current 
operator. We employ the form factor parametrization 
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decay mode 90% C.L. exp. bound on B la upper bound in units (mA/400 GeV)'* 

B+ ^ TT+r£+ 4.9 X 10-** \YebY:af + In^b^^r < 3-0 x 10"^ 

B+ -^ 7r+e±^T 1.7 x 10"" \YebY;/ + \Y^bY:af < 1.1 x 10"® 

B+ -^ K+e^fi"^ 9.1 X 10"* \YebY;,\'^ + \Y^bY:,f < 4.3 x lO"'' 

B+ -^ K+T^^j,"^ 7.7 X 10~^ \Y^bY;,\^ + \Yf,bY;,f < 5.7 x lO""* 

Table II: Limits on Y couplings coming from upper bounds on _B+ — > n{K)£^£'^ branching fractions, compiled by [23]. 

decay mode 90% C.L. exp. bound on Z? la upper bound in units (7tia/400 GeV)* 

8.0 X 10"* \YedY*a\^ < 1.9 X 10"^ 

1.1 X 10-^ \Y^,dY:f < 2.7 x 10"^ 
3.3 X 10-* \Y,dY*, - Y,,Y:f < 3.2 x IQ-^ 
4.0 X 10-* \Y^dY*s - Y^,Y:a\^ < 4.0 x 10"^ 

6.5 X 10-* 10.69 y^dy;^ - y^sY;^ < 1.3 x 10-^ 

Table III: Limits on Y couplings coming from upper bounds on r — > P£ branching fractions, determined at the _B-factories and 
compiled by [23]. 

for Pd = TT and Pg = K, where s = {p — p')^. The f^g form factors have been computed using QCD sum rules 
techniques and we employ the results of [39] . For the /^ q form factors we use a more recent calculation [40] . The 
B — > Kt^ 1.1^ differential decay rate can be written in a compact form by neglecting the small muon mass 
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(24) 
where A = A(l, mj^/rri^, s/m\) and A(a, 6, c) ^ a? + 1"^ + c^ — 2{ab + bc + ca) . For the modes without tau leptons in 
the final state one can neglect lepton masses completely, i.e. 

^^-(S ^ Ke^fi^) = Md!i^V^^3^A3/22^^(^)2 _ (25) 



ds (167r)'^?7i 



The modes with a pion in the final state can then be simply obtained from the above formula by replacing s with 
d and K with n . Integrating over the available phase space and comparing to the experimental upper bounds on 
B+ — >■ Tr{K)£~£'^ decays [23], we obtain the constraints listed in Tab. II. Finally we note that the corresponding 
rare K -^ 7r£+£'" decay modes are always less sensitive to the relevant Y entries compared to the rare leptonic 
^L,S^^+^'~ modes [21]. 

6. LFV semileptonic t decays 

These decays constitute important observables, uniquely sensitive to the third row of Y . Upper limits on their 
branching fractions have been set by the Belle and BaBar experiments. The width of the pionic channel reads 

where we have kept the leading powers of final state particle masses. Decay width for a channel with ^5 in the final 
state is obtained from (26) by replacing Y^dY*^ -^ YgdY*^ — YigY*^, /^ — >• \/2fK, and TO^r — >■ rriK- For the decay 
channel r — ^ /ijy we include amplitudes for both ss and dd components of 77 by replacing in Eq. (26) \YidY*^\ /^ — > 

\f^Y^dY*d + V^f^Y^sY*g\ , where /^''' are the decay constants of 77 through {dji^j^d + u7^75u)/-\/2 and 57^755 
operators, respectively. Following [41], we include the effects oirj — r]' mixing by using /^ = fq cos (j) and /^ = —fs sm<j) 
with phenomcnologically viable numerical values of fq = 1.07/^, fs = 1.34/^, and = 39.3°. With remaining 
numerical values /^ = 130.4 MeV [42], fx = 156 MeV [27], and the relevant 90% C.L. upper bounds on the branching 
fractions [23] we find a set of constraints shown in Tab. III. 



Nucleus V^^\mf'^'>] V^"\m^i:'^\ r.apturellO'^s- 



Au^i^ 



0.0396 
0.0974 



0.0468 
0.146 



2.59 
13.07 



Table IV: Data taken from Tables I and VIII of [43]. 



conversion in nuclei 



Four fermion effective Lagrangian (6) contains also the LFV terms {d'y^Piid){il^^PRe, ej'^Pufi}. Tlie most stringent 
bound on such interactions is expected from experimental searches for /^ — e conversion in nuclei. In order to derive 
the relevant constraints one needs to calculate the appropriate nuclear matrix elements of the above operators. A 
detailed analysis has been carried out in [43] . We can write the nuclear fi — e conversion rate as 



r, 



conversion 






l^(p) + 2y(")| 



(27) 



where the nuclear matrix elements y(P'"), calculated in [43] for titanium and gold nuclei are given in Tab. IV. 
Presently the most stringent bounds on S^e = rconvcrsion/rcapturo was set by the SINDRUM collaboration with 
bZ''^ < 4.3 X 10-12 [44] and S^t"^ < 7 x 10"" [45], both at 90% C.L. 
expressions we obtain the corresponding 1 a bounds 



Comparing these with our theoretical 



\Y,,Y;,\' < 1.9(20) X 10-13 {^^y froiii Au(Ti) . 



(28) 



Note that the same couplings also appear in the tt" — > e*/i^ decay branching fraction, whose expectation is thus 
pushed far below the current experimental upper bound of ^ IQ-^^. 

C. One-loop effects of A 

Next we turn our attention to observables which are affected by leptoquark couplings of A at the one-loop level. 
These are K — K and B — B mixing amplitudes, LFV neutral current processes like the radiative /x and t decays, 
as well as flavor diagonal observables, such as the anomalous magnetic moments of leptons or the decay width of the 
Z to bb pairs. With the exploratory nature of our study in mind, we do not consider nonlocal loop contributions 
due to the effective four-fermion Lagrangian (6), since such effects are constrained by the tree-level processes already 
considered in Sec. IIIB. The particular case of new absorptive contributions affecting Bg — B^ oscillations will be 
discussed in Sec. IV. 



1. eK and Atuk 




Figure 1: K — K mixing diagrams with leptons and A in the box loop. 
The SM result for the dispersive mixing matrix element, relevant for sk is [46] 



^^^12K — 



127r2 



fxTriMBx [XlmSoixc) + >H'n2So{xt) + 2XcXtri3So{xc,xt)] 



(29) 



\^k\ 


2.228(11) X 10"^ 


[23] 


AniK 


3.483(6) X 10-'^ GeV 


[2.3] 


0. 


43.5(7)° 


[23] 


fK 


0.1560(11) GeV 


[27] 


Bk 


0.725(26) 


[27] 


K, 


0.94(2) 


[48] 


ni 


l-31(ti) 


[49] 


m 


0.57(1) 


[46, 50] 


m 


0.496(47) 


[51] 



Table V: Experimental, nonperturbative, and perturbative parameters relevant for ex and ArriK observables. 

So is the Inami-Lim box loop function [47] and factors A^ = VisV*^ are the appropriate CKM weights. ExpUcit A^ 
contributions are eHminated using the Glashow-Ihopoulos-Maiani mechanism. Parameters rji, 772, and 773 account for 
the QCD renormahzation effects and are known to NLO (771,2) or NNLO (773) order. The decay constant fx and the 
reduced bag parameter Bk, both nonperturbative QCD parameters, are provided by lattice QCD calculations. Values 
of all the relevant experimental as well as theoretical parameters are compiled, together with their uncertainties, in 
Tab. V. The K — K transition is mediated also by box diagrams involving the A and leptons, as shown in Fig. 1, that 
generate an additional right-handed current operator in the effective Hamiltonian [52] 



n 



A 
AS=2 



1 



E^^ 



uYi, 



(c^fl7^Sfl)(dK7^s/i). 



(30) 



The dispersive mixing matrix element M12K induced by A is therefore 

1 



-'"l2_R' 



3847r2m^ 



JK^xBRm 



Yuy;, 



(31) 



where we have neglected the small QCD running effects from the A mass scale to the EW scale and simply use 
772 to describe the renormahzation group evolution of A contributions down to the hadronic scale. The observable 
measuring the CP-even component of the K^ mass-eigenstate, e^, is defined as the ratio of isospin singlet amplitudes 
of Ks(L) -^ TTTT decays 



eK 



A{Kl ^ (7r7r)/=o) 



A{Ks -^ {7nr)i=o) ' 
and is related to the imaginary part of the dispersive mixing amplitude as [48] 

6*"^= lmMi2K 



^K 



^/2 Attt-k 



(32) 



(33) 



Here /^mx is the measured mass difference between Ki^ and Kg eigenstates, while cj)^ is the superweak phase, given 
by (pe = arctan(2A77ii<: /AF/f ). The overall factor k,^ contains long distance corrections and uncertainties [48]. The 
resulting constraint on the Y couplings is then 



G|.77^^Im[A^77i5o(a;c)+At772S'o(a;t) + 2AcAt773S'o(a;c,Xf)] + 



m 



16771^ 



Re[^yMy;,]im[^y«y;,] 



12V27r2 A77ZK 



e^l = 1.57(7) X 10"'^ GeV 



(34) 



where on the right-hand side, we have combined the experimental and theoretical (hadronic) uncertainties by summing 
them in squares. Nonetheless, some theoretical uncertainty coming from the QCD renormahzation factors 771,2,3 still 
remains on the left-hand side. In the fit we allow them to freely vary within the intervals determined by their 
theoretical uncertainties (see Tab. V). 

The measured mass difference Arr^i^-, on the other hand, mostly probes the real part of the mixing amplitude 
M^i2K [53]. It receives potentially important contributions from SM long distance dynamics leading to large theoretical 
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uncertainties in its prediction [49] . Therefore we conservatively assume that the short distance contribution of M^j^ 
must be smaller than half the experimental value of Attik at 1 cr C.L.: 



^m-K^ReM-,^^^ 



f'^rriKBKmRe 



J2^mYI 



< 1.74 X 10"^'^ GeV. 



(35) 



The conservative assumption for the bound (35) allows us to neglect uncertainties of all the theoretical parameters 
and extract the following 1 a bound on the real part of the Y combination 



Re 



J2^idY;, 



< 1.1 X 10" 



rriA 



400 GcV 



(36) 



2. Bd — Bd and Bs — Bs mixing 



The time evolution of the B — B system is described by the average mass m, width F, and three mixing parameters 

|A^i2|, |ri2|, = -arg(Afi2/ri2). (37) 

All five parameters can be identified by diagonalizing the effective Hamiltonian 



i^cff = M 



-r = 




ri2 



(38) 



whose off-diagonal elements are defined as (-ffcff)i2 = {B I l-i^^~^ + nonlocal interactions I B) /(2m) and ri2 contains 
all on-shell contributions of intermediate poles. Heavy and light [H and L) mass-eigenstates are defined as 



\B 



H,L, 



= p\B)±q\B) , 



(39) 



and their eigenvalues are, in the appropriate limit |ri2| ^ |Afi2|, in turn connected to the measurements of Am and 
Ar as 



Am = ruH - niL = 2|Mi2| , 

Ar = rL-rH = 2|ri2|cos( 



(40a) 
(40b) 



If CP violation in decays is negligible, one can extract the phase (p from the semileptonic time-dependent CP asymmetry 



asi(i) 



T(B{t) -^ l+X) - T{B{t) -^ (rX) _ AF 
T{B{t) -^ l+X) + T{B{t) -^ e-X) ^ A^ 



tan( 



(41) 



The overall AF/A?tj factor renders this asymmetry very small. Measurements of a^j in the i?^ system have been 
performed at the B-factories and a world average [54] is consistent with zero, albeit with much larger errors than 
the SM predicted value. Direct measurement of a^j is not available, however D0 and CDF experiments [54] have 
measured the charge asymmetry of same-charge dimuon events coming from inclusive 6-decays, which is a linear 



M) 



M 



(.s) 



combination of a^j and a^j , thus allowing one to extract a^j . Especially the D0 measurements [13, 16, 17] point at 

an unexpectedly large mixing phase cj) and exclude the SM value of a^j with more than 3 a significance [18, 19, 55]. 
More effectively, one can extract the phase in the dispersive mixing amplitude M12 from the time-dependent CP 
asymmetry in decays of B {B) to CP eigenstates 



4p(t) 



TiBit)^f)-TiB{t)^f) 
F(i3(t)^/) + F(S(f)^/) 



?7/Im 



sin Am t , 



(42) 



where we have assumed a tree-level dominated decay mechanism, negligible CP violation in the decay, and also 
\p/q\ = 1. CP parity of the final state is denoted as ry/. To leading order in lFi2/Afi2|, F12 cancels out and one 
is sensitive to the phase of M12 through Im(p/q) = Im(Afi2)/lAfi2l. This phase is interpreted within the SM as an 
angle of the unitarity triangle sin 2/3 (sin2/3s) in the case of the Bd (Bg) system. Note that the weak phase of the 
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sin 2/3 
Arrid 

Arris 


0.673(23) 
0.507(5) ps-i 
17.77(12) ps"^ 


[54] 
[23] 
[23] 




0.275(15) GeV 
1.237(32) 


[27] 
[27] 


VB 


0.55(1) 


[46, 56] 



Table VI: Experimental, nonperturbative, and perturbative parameters relevant for Ba(^g^ — B^fgf mixing. 



-2/3(.) 



absorptive part is negligible in the SM and also difficult to enhance in most NP scenarios and thus (j)d(s) 
(we will comment on the recent study [52] of A contributions to the absorptive amplitude in Sec. IV). 

At this point we shall include as experimental constraints only the measurements of sin 2(3 and the mass splittings 
Am^j, Anis- We will address the allowed ranges of (ps in the fit part in Sec. IV. The SM prediction for the dispersive 
matrix element M12 in Bd mixing (with obvious replacements in the case of Bs mixing), is dominated by the short 
distance box diagrams involving the top quark 



^SM ^ ^k^flmBBB{V,,V;drVBSo{x,) 



(43) 



The three theoretical parameters here are again the perturbative QCD renormalization factor tjb, and the nonpertur- 
bative hadronic parameters fs and Bb- Box diagrams with A, analogous to the ones of -ftT — ^ in Fig. 1, can shift 
the value of Mf^ by 






3847r2r7r| 



fBrriBBBVB 



Y.^MY;t 



(44) 



where we again neglect the difference between the A and EW matching scales. Instead of using Am^ and Ato^ 
as individual fit constraints, we opt to trade Arrid for the ratio Arris /Am^; depending on the hadronic parameter 
^(= Bb^^/Jb^/Bbj^/Jb^) that can be determined reliably using lattice QCD techniques 



{VtbV*,)^Soixt) + (32mi G|rrr2^)-i (EiYe^Y;,)' 



iVt,V;^rSo{xt) + (32miG|m2^)-i (EeYuY, 



iVtbVt:)^So{xt) + (32miG|m2^)-i 



Y^^isY,] 



Am^ TUBg 2 
AiTid ruB, 

Arri 



22.5(12), 



67r2 



Gp-m'^r]BBB, /|^ ruB, 



3.53(39) X 10" 



(45a) 
(45b) 



On the right-hand sides, we have combined the experimental and theoretical errors in quadrature. The relevant 
numerical inputs are compiled in Tab. VI. 

On the other hand, in the sin 2/3 constraint all dependence on theoretical (in particular hadronic) parameters drops 
out 



Im 



So{xt){VtbV;,r + (32miG|m2^)-i (EeYuY, 



aJ 



So{xt)iVttV;,r + (32miG|,r«4)-i (E.YuY;,) 



= sin 2/3 . 



(46) 



3. Anomalous magnetic and electric dipole moments 

The electromagnetic interactions of an on-shell fermion can be parameterized in terms of parity conserving and 
parity violating form factors [57] 



F2 



r^ = Fi7^ + -^za^^'g. + F^a^'^q.-f^ + F^{2mq^ + 9^^)75 , 

2771^ 



(47a) 
(47b) 



12 

where q = p — p' . This is the most general form of the photon off-shell amplitude obeying the Ward identity of 
quantum electrodynamics 

g^^^ = . (48) 

Renormalized charge of a muon is — e and so Fi(0) = 1 exactly. A finite Fj,{Q) would signal a nonzero electric dipole 
moment in presence of CP violating phases in the renormalized vertex. i^4(0) is called the anapole moment. The 
form factor i^2, which is the source of the anomalous magnetic moment, enters in the gyromagnetic ratio as g — 
2(Fi(0) + ^2(0)). Comparing precise measurements of these form factors against theoretical higher-order predictions 
presents powerful tests of the SM and its extensions. In the recent years, the experimental result on the anomalous 
magnetic moment of the muon a^ = (g — 2)^/2 from BNL [10] has been about 3 a above theoretical predictions within 
the SM [11] 

a;=fP = 1.16592080(63) x 10"^ , (49a) 

a™ = 1.16591793(68) x 10"^ . (49b) 

Treating both experimental and theoretical uncertainties as Gaussian, we may identify the missing contribution to a^ 

6a^, ^ al^P - af^ = (2.87 ± 0.93) x 10"^ , (50) 

with the presence of NP. The leading A contributions to a^ with A and down quarks di running in the loop (Fig. 2) 



A ^-^-^ A rfj ^-'-^ rfj 

i ' d, \ i i f A \ £ 
^ i 4 • ^ ^ * 4 ^ 

Figure 2: Diagrams with A and down-quarks contributing to the lepton anomalous magnetic moments. 

are expected to be of the order ~ l/(47r)^ m? /m\ ely^il^ and have been previously computed in [58]. We reproduce 
the magnitude of a^ of [58] , however with an opposite overall sign 

" ^ Y^ y^ \Yt,i\'^[QAfA{xi) + Qdfd{xi)], Xi^^ml Jm^. (51) 



^ i=d,s,b 

Here the charges are QA.d = 4/3, —1/3 while fA,d are the loop functions 



^ , , 2x3 + 3x2 -6x2 log a: - 6a; -fl 

fA{x) = ^^^-^p , (52a) 

-x^ +6x2 - 6xlogx - 3x - 2 
M-) 6(^^3T)i • (52b) 



In the limit Xi — )■ the result becomes 



^ i=d,sM ^ i=d,s,b 

If we now saturate ^a^ with a^ we find that a non-zero magnitude is preferred for a combination of the second row 
elements of Y 

Y: \Y,.\' - (4.53 ± 1.47) X 10'^ X ^ = (6.45 ± 2.09) x "f . (54) 

i=d,s,b A" V J 

We will further explore the possible correlations of such effects with other constraints in Sec. IV. 



13 



decay mode 90% C.L. exp. bound on B la upper bound in units (mA/400 GeV)* 



At-)- 67 


2.4 X 10"^^ 


r-> ^7 


4.4 X 10"* 


r — > e-y 


3.3 X 10-* 



YeiY;,\^ < 4.6 X 10-* 

2 ^ /I o s^ 1 n — 3 



IE.=d,.,,i;.«n*r< 4.8x10 
IE,=,,.y«n*r<3.6xlO-« 



Table VII: Limits on Y couplings coming from upper bounds of LFV radiative lepton decay branching fractions, taken from [23, 
59]. 

On the other hand, applying expression (54) to the electron case and requiring that af be smaller than the 
experimental uncertainty, we find a 1 ct bound on the first row of Y 



E 

i—d^s,b 



Y 



< 



X 10" 



= 54x 



TO/ 



(400 GeV)2 



(55) 



where we have used the experimental uncertainty estimate of u^j^'^p = 2.8 x 10" [23]. 

Finally, we note that due to the Hermitian structure of Y contributions to the EM interactions of quarks and 
leptons, no electric (or chromoelectric) dipole moments of either quarks or leptons are generated at the one loop 
level, regardless of the phases present in Y . Furthermore, even at the two loop level, non-zero contributions can only 
originate from mixed W — A loops. However, since A interactions are purely right-handed, such contributions are 
necessarily suppressed both by CKM factors and by insertions of the light quark or lepton masses. Therefore we do 
not consider them further. 



4-. Flavor violating radiative decays 



The computation of A contributions to the LFV radiative muon decay is analogous to the magnetic moment 
diagrams in Fig. 2 and results in the effective Lagrangian 



C 



M^e7 647^2^2^ 



E 

i—d.s,b 



V Y"^ 

^ PI J- 1,. 



e{a^''F^,){m^,PL + m,PR)fi. 



(56) 



The decay width of /i — > ej is then given by 



40967r4m^ 



E 

i^d,s,b 



Y F* 



(57) 



The above expression can also be applied to the LFV decays of the r, with obvious replacements in Y indices and 
masses. Inequalities following from upper limits on branching fractions of ^ ^ £'7 are shown in Tab. VII. Consequently, 
measurements of these flagship LFV processes impose strict requirements on the structure of Y ^ namely they require 
that rows of Y are approximately orthogonal. 

On the other hand, the analogous constraints coming from the quark sector radiative decays are much weaker. 
The prominent example of & — )■ 57 has recently been analyzed in [55], where it was found that this decay is not very 
sensitive to the relevant A interactions, which contribute at one loop through the insertion of the Lagrangian (6). In 
particular, the Anig constraint (45b) yields consistently stronger bounds on the same combination of Y elements for 
the experimentally allowed range of A masses. 



5. Decay width of Z ^ bb 



The experiments running on the LEP 1 collider performed precise measurements of the relative widths oi Z ^f hh 
and Z — >■ hadrons. In particular the experimental value of 



Rh 



T{Z -^ bb) 
T{Z -^ hadrons) ' 



(58) 
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b \ b 

Figure 3: Diagrams with A and leptons modifying the Zbb vertex. 



is according to the PDG [23] in good agreement with SM predicted value 

i?°''P = 0.21629(66) , 
Rf^ = 0.21578(5) . 

The SM tree-level amplitude for the Zbb vertex is 



A' 



= ^gz 



9r-Ar + glAL 



al 



SIT? Ow , gl 



■Al(r) = Z^b^f,Pnj^:jb, 
-l/2+-sm''0w, 



where gz — g/ cos 6w- New contributions of A change 5^ to qr ~ gR + ^gR, where 



^gt 



■2n EAYibl' 1 
sm uw 
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(59a) 
(59b) 

(60) 



xl - 2xz + \og{xz){3xl - 6xz + 61og(l + xz)) 



(47r)2 6a;| 
1 + 4/2(2x2 - 4) + 6U2{-xz) - 1^(34 -6xz + 6 log(l + Xz) 



(61) 



/m^. Calculational details along with functions /i and /2 are given in the Appendix A. Taking into 



and Xz — '"'z/ '"'A 

account higher order SM corrections, the relative shift in Rt due to such NP contributions can be written as [60] 

SRb = 2R^^(l - i?SMv9i Re((55L) + 9% Re{6gB.) 



{9l? + {gl? 



(62) 



One can check that the shift Re(J(7/j) given in Eq. (61) is negative and Sql = and consequently any contributions 
of A necessarily worsen the agreement between theory and experiment. If the discrepancy should be smaller than 1 a 
the following constraint has to be met 



^|ya|2<5.60( 



TOA 



400 GeV/ 



6.73 



mA 



.400 GeV'y 



2.02. 



(63) 



In derivation of this bound we have approximated Sg^ by a polynomial in variable mA and employed sin 9w ~ 
0.231 [23]. 

On the other hand, the forward-backward asymmetry in bb production as measured at LEP exhibits a 2.1 a tension 
with the SM EW fit. Since A contributions only affect the right-handed effective Zbb coupling {5gR) they cannot 
fully reconcile this tension [60] and we do not include this observable in the fit. 



IV. GLOBAL FIT OF THE LEPTOQUARK COUPLINGS 

In this section we perform a global fit of Y to all the observables listed in Sec. Ill, while we keep fixed ttia = 400 GeV. 
We resort to a x^(^) statistic that we minimize to find the point Xmin = X^(^bcst), where Y is by definition in best 
agreement with all the constraints, x^(^) is written as a sum of Gaussian contributions of observables Oi 



\y) = j: 



o^ 



o. 



prediction 



{Y) 



^cff\ 



(64) 
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Values of O^^^ and a^^ are central values and errors, read-off from right-hand sides of constraining equations in the 
preceding sections, whereas Of^" "^ '°" are the corresponding predictions in terms of Yij, i.e., the left-hand sides of 
constraints, in the language of Sec. III. Majority of 0°'^^ are upper bounds which are modeled with x^ centered at 
zero. This is achieved in Eq. (64) by setting 0°'^^ — and a^^ to the derived 68% C.L. upper bound. Although not 
explicitly shown here, the x^ function (64) depends also on the 4 Wolfenstein parameters of the CKM matrix (they 
are present in meson mixing constraints) and we treat them on the same footing as Y. We add to (64) a Gaussian 
chi-square term which guides the CKM parameters to follow probability distribution of Eqs. (4) and (5). 

Statistical interpretation of the value of x'^^ i-6- the goodness of fit, is performed using the standard x^ probability 
distribution with appropriate number of degrees of freedom, A'dof = -^obsorvabics — -^parameters = 36 — 21. To find the 
allowed range of a single matrix element \Yij\, its phase, or a function of several Y elements, denoted in the following 
generically as z{Y), we minimize x^ with z{Y) fixed to some chosen value zq. Then all values Zq, where 

Tam[xHY).(Y)=z„] - Xmin < 1 (4) (65) 

form the 68.3 (95.45) % C.L. interval for the parameter z. To find confidence level regions in two-dimensional scans 
(with two fixed quantities z{Y), w{Y)) we utilize the ^^-distribution with 2 degrees of freedom. The difference 
^X^{zo,Wo) = min[x2(y)=2o,„(i')=u.o] ~ ^min ^^^' pomts {zo,wo) in the Na C.L. region in this case is 

1 - exp [-Ax^{zo,wo)/2] < erf (^^ . (66) 

A. Structure of Y 

In a trivial case, when we set Y — to recover the SM, we find Xinin — 12.5 = 9.5a + 1-5ckm + Q-^Am^ + ■ ■ ■ with 
a dominant contribution from the a^ anomaly. If we let Y take any value we find a global minimum Xmin = 2.5 = 
l-8cKM + 0-4:Am, + ■ • • for 15 degrees of freedom, which signals a very good agreement of all predictions with the 
considered observables. In particular, the best point perfectly resolves the anomalous magnetic moment constraint a^ 
and slightly improves quark fiavor constraints. The allowed 1 and 2 a ranges of Y matrix elements are shown below 

(< 1.4x10"^ < 8.7x10^5 < 4.2x10"^ \ 

< 3.6 X 10-3 U [2.1, 2.9] < 3.6 X 10-3 u [2.1, 2.9] < 6.2 x IQ-'* U [2.3,2.7] , (67a) 

< 5.6 X 10-3 < 8.1 X 10-3 < 9.6 x 10-3 J 

I < 2.2 X 10-6 < 1.4 X 10-"* < 6.6 x 10-* \ 

\Y^'^''^\ e < 5.6 X 10-3 u [1.5, 3.3] < 5.6 x 10-3 u [1.5,3.3] < 9.7 x IQ-* U [1.6, 3.2] . (67b) 

^^ < 8.9 X 10-3 < 1.4 X 10"2 < 1.5 X 10-2 J 

Couplings to the electron are strongly suppressed, while couplings to the muon (the second row of Y) can take values 
of order 1, in order to satisfy the a^ constraint. In the last row, elements Yts and l^b can also be of order 0.01 at 1 cr 
C.L. We find some interesting correlations between the second and third row elements, shown in Figs. 4 and 5. 

We find three distinct regimes in the second and third row (Figs. 4 and 5), depending on which element in the 
second row is large. Pictorially, these hierarchies are possible 

^0 o\ /o o\ /o o\ 

■ 00, OBO, OOB. (68) 

V' • V V' • V \' • •/ 

Here ■ stands for order 1 element, • for (at most) order 0.01 element, while we neglect elements which are < 10-3. 
This particular hierarchy is enforced by a nontrivial a^ constraint that requires at least one large element in the 
second row, while stringent upper bounds from LFV processes exclude the possibility of having two elements of order 
1. 

We also identify the observables, which are most constraining for each element in Y. We do this by registering 
the maximum increase in each individual observable contribution to x^ when a single Y element is changed from its 
best-fit value. In the first row, all the most stringent constraints actually also involve the y^g elements: Kl —>■ /x^e^ 
{Ygii), /i — ^ 67 (Yes) and B — ^ ■nji^e^ (i^eb)- In the near future, we can expect some significant improvement at least 
for Yes from the MEG experiment [61, 62]. The best-fit regions around 0(1) values for the second row elements are 
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J3 



0.02 



0.002 



O.OOO: 




3002 0.002 0.02 0.2 



'^sl 



0.002 




0.02 



1.002 0.02 0.2 

|F„| 



Figure 4: Correlation between elements in the second (IVpsI and IV^tl) and third {\Yts\ and |yr6l) row. Dark green region is 
the 1 a contour, while the two lighter green contours are 2 and 3 a, respectively. 



o.oo: 




0.02 



0.0002 0.002 0.02 0.2 



Figure 5: Correlation between diagonal Y elements in the muon and tau rows. Dark green region is the 1 a contour, while the 
two lighter ones are 2 and 3 a, respectively. 



mostly determined by the observed discrepancy in the a^. Other relevant observables, that constrain their values 
in the <C 1 regions are: AniK and Ks — >■ /x+zz" (F^d, Yfis), and B — >• Xsi~^i~ {Yf_ib)- Unfortunately, due to the 
theoretical uncertainties which dominate the precision of the first two observables, a significant improvement in the 
foreseeable future can only be expected for the constraint on y^if, from the Super Flavor factories (SFFs) [63, 64]. 
Finally, the constraints on the third row of Y are dominated by LFV tau and B decays: r — )> /iTTo and t ^ firj (Yrd), 
B — > Kr^fi^ (Yrs), and B -^ Kr^fi^ and t — > ^Ks {Yrb) ■ Again SFFs are expected to yield improved bounds on 
these LFV observables. Finally we note the fact that bounds on most of the elements of Y are dominated by rare 
decays and the a^ which all exhibit a similar scaling dependence on the A parameters {Y/m/\). This points towards 
an approximate linear scahng of the fitted Y element values with the A mass and allows for simple reinterpretation 
of the derived limits at A masses away from the reference value ttia = 400 GeV. 
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B. Comment on tension between B{B ^f tv) and sin 2/3 

We can redo the tree-level CKM fit, described in Sec. Ill A, replacing the \Vub\ value from Eq. (3) with a constraint 
coming from a world average of Belle and BaBar measurements [65-68] of B{B^ — > r+i^) = (1.68 ± 0.31) x 10^"* [19]. 
The observable cannot be directly affected by A contributions and is given in the SM by 



GpTng+ml 

Sn V 'm 



g(J3+^rV)= "^ :r M l-^) lK.|V|r^+. (69) 



The main theoretical uncertainty due to the lattice QCD estimate of the relevant decay constant /b = 193 ± 
10 MeV [27] is at present subleading compared to the experimental error, but we nevertheless combine them in 
quadrature. The best-fitted values of the CKM parameters are then 

A = 0.22538(65) , (70) 

A = 0.799(26) , 
p = 0.162(90) , 
77 = 0.528(64) . 

The quality of the fit of the CKM from tree-level observables is exactly the same as in Sec. Ill A, however the central 
values of p and especially 77 are significantly higher than before. This is expected since the tree level fit of p and 
7] parameters is not over-constrained. We can also repeat the global fit of couplings Y, this time with a Gaussian 
chi-square term for the CKM matrix corresponding to Eq. (70) and the underlying correlation matrix. The best fit 
point with x^ = 9.5 relaxes the tension in the CKM by changing the Bd — Bd phase (sin 2/3) but at the price of not 
resolving the a^ anomaly at all. Another, slightly shallower, minimum with x^ = 10.6 achieves just the reverse - a^ 
is perfectly satisfied while the tension in the CKM persists. 

A qualitative explanation goes as follows: a new phase in Bd — Bd mixing can be generated by either Y^dY*ij, Yfj.dY*f, 
or Y-rdY*^. Large YedY*f^ and Yf^dY*^, are ruled out by the strong B -> irl^i^ constraint. Thus for large enough YrdY*fj 
either (i) Y^d or (ii) Y-rb is at least ~ 0.1. In turn we form uncomfortably large products of (i) IVd^^g or (ii) YrtY^^q 
for q ~ d, s, b, where Yf^q{^ 1) is large to explain a^. First possibility {q = d) is incompatible with (i) r — >■ Tr^/i or (ii) 
Bs — > Tp branching ratios, second one with (i) t — > K° p or (ii) B — > Krp, and the last one with (i) Bd — > rp or (ii) 

T — > p^. 

Such worsening of the overall agreement of observables with the model can already be anticipated from Eqs. (67) 
which clearly state that large contributions to Bd — Bd mixing are disfavoured. On the other hand, if one ignores the 
a^ constraint, then this model can sufficiently affect the phase of the Bd — Bd mixing amplitude to be consistent with 
a large B{B — > tv). 

C. Comment on CPV in the Bs system 

Next, we address the question whether contributions of A can enhance the phase of the dispersive amplitude 
in the Bs — Bs system, or even modify the absorptive part by on-shell charged leptons in box diagrams (Fig. 1). 
According to [52], which studied tt absorbtive contributon to mixing amplitude, one should have yj-^y^j ~ 0.1 for 
?7iA = 250 GeV. However, the set of observables we have included in the analysis of Y couplings (in particular a^) 
forces the tt and pp states to couple very weakly to Bs 

|i;,r;b I < 0.0015 (0.0021), (71a) 

\YrsY;^\ < 0.9 X 10"'' (4.1 X 10"'') . (71b) 

Bounds are 1 a (2 a) C.L. (see also Fig. 4) . While [i^^s^utl i^ directly constrained by the B — > Xgi'^i^ rate, there is 
presently no direct bound on the magnitude of Yt-sY*jj, so the constraint (71b) is directly linked to the explanation 
of the a^ anomaly. Using the above 2 a upper bounds in Eq. (44) we find that the dispersive A amplitude with tau 
(muon) in the box is five (four) orders of magnitude smaller than the SM contribution. 

The approximate scaling of the most important constraints with Y/m^ provides a robust bound on the absorptive 
contributions to the neutral meson mixing amplitudes, excluding any significant modification of AF^, provided we 
require the resolution of the a^ anomaly. We also note in passing that the smallness of new absorptive NP contributions 
is required in general by the measurements of the lifetime ratios of B mesons, semileptonic branching fractions, and 
the average number of charm quarks in B decays (a recent discussion can be found in [G9], see also [19] and references 
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Figure 6: Correlation between the absolute value of YtsY*i, and its phase {ATg{YTaY*i,)) in absence of the a^ constraint. Dark 
green region is the la contour, while the two lighter ones are 2 and 3a, respectively. 

therein) . On the other hand, the maximum allowed relative A contributions to the dispersive parts scale quadratically 
with rriA ■ In this way dispersive A amplitudes comparable in size to SM contributions in Bs — Bs mixing observables 
are only reached at masses well above 1 TeV, where the relevant Y couplings are no longer perturbative. Thus we 
find no possibility to simultaneously affect a^ and the Bs system observables with A contributions. 

In absence of the a^ constraint, the bounds on YTgY*^^ are significantly relaxed and are dominated by AttIs and 
Ams/Arrid (see Fig. 6). |lVsXr*{,| values of order 0.1 are allowed, however at the expense of fine-tuning the phase 
ATg{YrsY*f^) in order to obtain the right destructive interference with the SM contributions to the i?s,d mass differences. 

V. GUT IMPLICATIONS 

A. Framework 

The color triplet leptoquark (3, 1,4/3) emerges naturally in a theoretically well- motivated class of grand unified 
models. We will first demonstrate this in a framework of the SU{5) gauge group — the simplest group to encompass 
the SM gauge symmetry — and then proceed to discuss how and where it appears in the 5*0(10) setup. 



1. SU{5) setup 

The matter of the SM is assigned to the 10- and 5-dimensional SU{5) representations, i.e., 10^ = (1, 1, 1) 
(3,1,-2/3) © (3,2,1/6) and 5, = (1,2,-1/2) ® (3,1,1/3), where i{= 1,2,3) denotes generation index [70]. This 
assignment dictates that the charged fermion masses and the entries of the CKM matrix originate, at the tree-level, 
through the couplings of the matter fields to the 5- and 45-dimensional Higgs representations only [71]. It has 
actually been shown that the phenomenological considerations require presence of both [72-77] . It turns out that the 
color triplet leptoquark is a part of the 45-dimensional representation. Namely, the relevant SM decomposition reads 
45 = (Ai,A2, A3, A4, As, Ag, At) = (8,2,1/2)0(6, l,-l/3)®(3, 3,-1/3)0(3, 2, -7/6)®(3,l,-l/3)e(3,l, 4/3)® 
(1, 2, 1/2). The color triplet thus appears in any SU{5) framework that relies purely on the scalar representations for 
the charged fermion mass generation. 

Relevant contractions of the 45- and 5-dimensional Higgs representations, i.e., 45 and 5, with the matter fields, 
are (Yi)yT0i5j45*, {Y2)ijl0il0j45, {Y3)ijl0^5j5* and (F4)ijT0,jl0j5, where Ya, a — 1,2,3,4, represent arbitrary 
Yukawa coupling matrices in fiavor space. The charged fermion mass matrices at the unification scale accordingly 
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read 

Md - -Y,vl,^]^Y:,vl (72) 

Me = iY^vl^-^Y^vl (73) 

Mu = 2V2{Y2-Y^)vi5-^{Yi + Yl)vz, (74) 

where (5^) — v^jy/l and (45p) = (45|^) = (453^) = ^45/1/2 represent appropriate vacuum expectation values. Note 
that 5 = 5", 45 = 45^^ and |w5lV2 + 12|w45p = f ^ where a, /3, 7 = 1 ,.., 5 represent SU{h) indices and w(= 246 GeV) 
stands for the electroweak vacuum expectation value (VEV) . (The VEV result has been introduced for the first time in 
Ref. [G] and corrects the normalization presented in Refs. [7G, 78].) In SU{5) there could be an additional contribution 
to V from an SU{2) triplet scalar [79] but that contribution is supposed to be suppressed by a large symmetry breaking 
scale [80] and we accordingly neglect it. We also assume that both U5 and V45 are real for simplicity. 

In order to have consistent notation we identify Ag with A in what follows. The lepton and baryon number violating 
Yukawa couplings of the triplet A to matter in the fermion mass eigenstate basis in the SU{5) framework are already 
given in Eq. (2) if one makes the following identifications: Y = _e]^Yi£)|j and g = 2\fw\\Y2 — ^2^]C^i?/ Here, Er^ 
Z?fl and Vji represent appropriate unitary transformations of the right-handed charged leptons, down-quarks and 
up-quarks. Our phenomenological study primarily relates to Yukawa couplings of A to the down-quark and charged 
lepton sectors. Clearly, these low-energy constraints on the leptoquark couplings to the matter could allow us to place 
constraints on the very Yukawa couplings and associated unitary transformations that show up in the charged fermion 
mass relations. These, on the other hand, might be pivotal in addressing the issue of matter stability [81]. 

Note that the antisymmetric nature of the color triplet couplings to the up-quark sector in Eq. (2) is dictated by 
the group theory and is not affected by any change of basis. In other words, any unitary redefinition of fermion fields 
would preserve this property. We insist on this point for the following two reasons. Firstly, this is important since 
it is this unique feature of the A couplings to the up-quark sector that is responsible for an absence of the leading 
contributions towards proton decay due to A exchange [78] . Secondly, if, for some reason, Yi is a symmetric matrix, 
there would not be any coupling between A and the up sector. In other words, all gij elements in Eq. (2) would be 
zero. If that was the case, A would not mediate proton decay. In fact that can happen, for example, if the scalar 
leptoquark A originates from an 50(10) setup as we discuss next. 

Z. 50(10) setuv 

Recall, one generation of the SM matter in the 5*0(10) framework is embedded in a single 16-dimensional repre- 
sentation. The allowed contractions of the matter fields to the Higgs sector, at the tree-level, are (Yio)ij 16^16^10, 
(Yi2o)ijl6il6jl20 and (Yi26)ijl6il6jl26, where 10, 120 and 126 are the scalar representations that all contain 
states with the quantum numbers of the SM doublet [71]. Here, Yio(= ^10) : ^120 (= ~^i2o) ^"^^ ^126 (= ^126) ^^P" 
resent complex Yukawa coupling matrices. As it turns out, the 45-dimensional representation of SU{5) is found in 
both the 120- and 126-dimensional representations [71]. The former one couples antisymmetrically to matter, thus 
preserving the absence of the leading contributions towards proton decay due to A exchange [78]. The latter one, 
on the other hand, couples symmetrically to matter. So, if A originates from the 126-dimensional representation of 
5*0(10), it will not couple to the up-quark sector at all. Consequently, there will be no proton decay signatures related 
to A exchange in that case. Again, these properties are dictated by gauge symmetry and are preserved regardless of 
any redefinitions of the charged fermion fields. (Note that our findings on the absence of the up-quark sector couplings 
do not agree with the conclusions put forth in Ref. [9] for the 50(10) case and in Ref. [82] for the SU{5) case.) 

The relevant mass matrices for the down-quarks and charged leptons in the 50(10) framework are 

Md = -Yi26V*u6-^^wvlo + Yi2o{v[*2o + <;o), (75) 

Me = 3Yi26Wi26- 2^10^10 + ^120 M*2o-3<2*o)> (76) 

where wio, wi26j ^^^20 and i'"20 represent VEVs of the doublet components of the appropriate scalar representations. We 
will assume that the VEVs are real when needed for simplicity. (See Ref. [83] for exact normalization with respect to 
the SM VEV.) Clearly, the observed mismatch between the charged lepton and down-quark masses requires a presence 
of either 120 or 126, or both representations in the case without the 10. The color triplet hence must appear in any 
50(10) framework that relies purely on the scalar representations for the charged fermion mass generation. 
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We opt to start our analysis within a particular class of SU{5) models having in mind that the same procedure 
can be carried over into an SO{10) framework with appropriate modifications. In fact, towards the end of the next 
section we also address the 150(10) setup viability in view of its compatibility with phenomenological constraints on 
the couplings of the light colored scalar to the matter fields. 

B. Numerical Analysis 

Our goal is to consistently implement all available constraints on the color triplet couplings to the down-quarks 
and charged leptons in order to study implications for the charged fermion Yukawa sector within a particular class 
of grand unified models. These models rely solely on the scalar representations in order to generate charged fermion 
masses. 

We first single out a simple SU{5) setup with the 5-, 24- and 45-dimensional representations in the Higgs sector [72- 
77] and one 24-dimensional fermionic representation [77] to generate neutrino masses via combination of type I [84- 
88] and type III [89, 90] seesaw mechanisms for definiteness. We resort to this model since it has been explicitly 
demonstrated that it predicts proton decay signatures that are very close to the present experimental limits on the 
partial proton decay lifetimes for the mass of A in the range accessible in collider experiments [2]. (The model is a 
renormalizable version of the scenario first proposed in [91] and further analyzed in [92, 93].) Moreover, it shares the 
same mass relations given in Eqs. (72), (73) and (74) with all other SU{5) scenarios that rely on the use of the 5- 
and 45-dimensional scalar representations. 

We start with the following relations that are valid at the unification scale 

eIDlM^''^^ = {-^eIY3D*j,v,~Yv^5), (77) 

Mf'-'^ElD*^ = {-]^ElY^D*j,v^ + iYvi^), (78) 

where Mj^^^ and M^^^ are diagonal mass matrices for down quarks and charged leptons, respectively. Our con- 
vention is such that Me, — DlMj^^^D^ and Me = E^M^^'^'^E]^, where D^ and E^ represent appropriate unitary 
transformations of the left-handed down quarks and charged leptons. Note that our phenomenological considerations 
yield constrains on the form of Y that are valid at low energies only. It is thus essential to propagate constraints on 
{Y)ij, i,j = 1, 2, 3, as well as the entries of Mj^^^ and M^^^ to the GUT scale to extract accurate information on W45 
and unitary matrices Ej^Dj^ and Ej^D'^. 

Again, the phenomenological bounds we derive constrain the matrix Y appearing on the right-hand side of a relation 

eIjDlM'^'"^ - m'^'^^eId*^ = -4ru45. (79) 

What is not known are the overall scale of the right-hand side set by W45 and the unitary transformations given by 
Ej^Dl and E'[D*j^. In order to perform numerical analysis and implement inferred bounds we first parametrize EJ^Dl 
and EJ^D*j^ using a generic form 
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where Sab = sin^a;,, Cab = cosOab, and ai, i = 1, ..,6, are phases. We then randomly generate the total of nineteen 
parameters and check whether the left-hand side of Eq. (79) satisfies all phenomenological constraints. (We also vary 
the four parameters of the CKM matrix as well as rji, 772 and 773 — QCD parameters entering K — K mixing — in 
order to have consistent constraints on the Y entries as described in Sec. III.) This process is repeated until the 
available parameter space is thought to be exhausted. The down-quark and charged lepton masses at the GUT scale 
are considered as input and the relevant values we generate and use within this particular framework are given in 
Table VIII. The GUT scale is taken to be Mqut = 10^® GeV and we only propagate and use the central values for 
the down-quark and charged lepton masses. 

Note that the need to accommodate experimental results on a^ basically sets the scale for the Y entries. To be 
precise, it requires that X]i=i 2 3 l^2i|^ satisfies Eq. (54). This in turn should fix the value or range of allowed values 
of Viz since the scale of the left-hand side of Eq. (79) is set by the known fermion masses. One can then use this 
information to determine v^ via |d5|^/2 -|- 12|w45|^ = ^^- To be conservative we not only vary X)i=i 2 3 l^2i|^ within 
the 1 a and 2 a ranges but accommodate for the effect of the RGE running of our constraints from the low scale to 
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running mass at Mz 


running mass 


at Mgut 


mt{Mz) 


= 2.89 ± 0.11 GeV 


mi,(MGUT) = 


0.782 GeV 


ms{Mz) 


= 56±16MeV 


ms(MGUT) = 


19 MeV 


md{Mz) 


= 3.0±1.2MeV 


md(MGUT) = 


1.1 MeV 


m.r{Mz) 


= 1746.45l°;^8 MeV 


^-^(Mgut) = 


1561.4 MeV 


m^{Mz) 


= 102.72899(44) MeV 


^.^(Mgut) = 


= 9 1.84 MeV 


m^{Mz) 


= 0.4866613(36) MeV 


me(MGUT) = 


0.4350 MeV 



Table VIII: Input parameters for the relevant fermion masses at the Mz scale and the corresponding values at the GUT scale 
(A^GUT = lO^'' GeV) in a non-supersymmetric framework. 

the grand unified scale. We take that effect to be within the bounds set by the following scaling factors: 1.1-3.7. 
These scaling factors correspond to the maximal changes in the charged lepton and down quark masses as they are 
propagated from low scale to the GUT scale. Again, we take the GUT scale to be Mqut = 10^^ GeV for simplicity. 
(The exact dependence of the GUT scale on the scalar particle mass spectrum within this particular SU{5) model is 
known and has been worked out in detail in Ref. [2]. The change in the GUT scale or, correspondingly, the scalar 
particle mass spectrum also affects propagation of fermion masses but that effect is rather small for the scenario when 
A is light as the GUT scale is then limited within a very narrow range [2].) 

The upper limit on U45 which we obtain by randomly choosing the entries of Ej^D^ and Ej^D'^ is shown in Fig. 7. 
Clearly, the bound should drop as mA grows since V45 needs to compensate the growth of the appropriate values of 
Y that satisfy the a^ constraint of Eq. (54). For practical purposes, we generate this conservative limit when only 
Eq. (54) is satisfied for a finite set of fixed values of niA- These correspond to dots in Fig. 7. In our numerical study 
we limit the toa range due to the existence of both the lower and upper bounds on its value. The lower experimental 
bound on ttia comes from direct experimental searches. The most stringent one originates from dedicated searches for 
pair production of leptoquarks in pp collisions at LHG and it reads ttia > 384 GeV [94] (mA > 422 GeV [95]) for the 
so-called first-generation (second-generation) leptoquarks assuming these decay exclusively to an electron (muon) and 
a hadronic jet. While these bounds are not necessarily applicable to our framework, since A can also decay to a top 
quark and a hadronic jet, we have verified that the corresponding branching ratio is always below 30% in the region 
of parameter space where A resolves both the tt FBA and the a^ puzzles. The upper bound on ttt-Aj on the other 
hand, originates from perturbativity constraints on entries of Y that should not exceed v47r- We find that bound to 
be ruA < 560 GcV. 

After an extensive numerical study we fail to generate a single satisfactory solution to all the constraints using 
Eq. (79) as a starting point. We trace the difficulty of finding a viable numerical solution to the facts that (i) the 
down-quark and charged lepton sectors do not exhibit a strong mass hierarchy that is present in the up-quark sector 
and (ii) the misalignment between the masses of the down-quarks and charged leptons that belong to the same 
generation is sufficiently large to prevent necessary cancellations. For example, a generic form of the left-hand side in 
Eq. (79) can be represented as follows 

/o o\ /go •\ 

(81) 

Here ■ (•) stands for an order rrir (mb) element. Clearly, the only potentially viable scenario for this form to 
describe matrix Y, pictorially given in Eq. (68), would be the one where the 23 element dominates. The 31 and 32 
elements should accordingly be suppressed by effectively setting the angles ^13 and 623 from Ej^D*^ to zero. This, 
however, leaves the 33 element on the left-hand side of Eq. (79) to be proportional to ra-r — mh(iJ|jZ?i)33. As 
77iT-(AfGUT) ~ 2TO;,(AfQUT) in the scenario at hand and |(i?]jDi)33| < 1, the absolute value of the 33 element turns 
out to always be greater than the absolute value of the 23 element, in contrast to what is needed. One could try to 
see if there is a possibility to have a satisfactory numerical solution within the supersymmetric framework where, for 
example, the mismatch between h and r varies a lot with the change in the tan/3 parameter. This scenario, although 
it does help in suppressing the 33 element, also fails due to the difficulty to accommodate small enough elements in 
the 1-2 block of the left-hand side of Eq. (79). Namely, once the freedom to set the 13 and 33 elements to be small 
by tuning the angles in EJ^Dl is used there is not enough parameters left over to tune the 1-2 block to the desired 
form. For example, since me(MGUT)/wf,(MGUT) ^ "^(j(MGUT)/"i&(MGUT) ~ 10~^, the 11 element is always bigger 
than the required limit of 10~^. In short, the SU{5) scenarios with a light triplet scalar that rely on the use of the 5- 
and 45-diniensional scalar representations to generate charged fermion masses at the tree level fail to accommodate 
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Figure 7: Upper bound on ^45 as a function of ttia- Data are generated for a discrete set of ttia values that are shown as dots. 
The curve is an interpolation that carries an in^^ dependence. 

the Yukawa structure needed to explain the a^ puzzle while satisfying all other phenomenological constraints. 

Let us now discuss implications of our findings with respect to their compatibility with the most commonly encoun- 
tered 50(10) scenarios. Recall, the only representations of SO{10) that could, at the tree-level, yield contributions 
to the charged fermion masses are the 10, 120 and 126. And, as we have pointed out in Sec. VA, A can originate 
from either 120- or 126-dimensional representation of 50(10). 

If A is part of the 126-dimensional Higgs it would not couple to the up-quark sector since the relevant couplings to 
matter are symmetric whereas A needs to couple in an antisymmetric manner to the up-quarks. If, in addition to the 
126, one uses a 10-dimensional scalar representation to generate the charged fermion masses the corresponding mass 
matrices will all be symmetric. This, on the other hand, changes the transformations that generate mass eigenstate 
basis from bi-unitary into congruent form. This significantly reduces the number of free parameters yielding the 
following mass relation 

C/M^'^s _ M^iag^* ^ -AYV126, (82) 

where U = E'j^Dn, Eji = E^ and Dj^ — D^. This relation also corresponds to the SU{5) scenario when all Yukawa 
matrices in the down-quark and charged lepton sectors are symmetric. Obviously, this case is much more restrictive 
since we have only one unitary matrix U to vary. It is thus clear that this scenario cannot be viable if we implement 
all the constraints on the form of Y. Hence, the case when Yukawa couplings in the charged lepton and down-quark 
sectors are symmetric, including the case with the 10- and 126-dimensional scalar representations in 50(10), is not 
compatible with possibility to have light A as an explanation for observed anomalies. 

The scenario with the 10- and 120-dimensional representations in the Higgs sector is also not realistic. In fact, that 
scenario resembles the SU{5) scenario that proved to be inadequate to accommodate the form of Y matrix. Moreover, 
the 10- and 120-dimensional representation scenario cannot explain observed fermion masses as was demonstrated in 
the low-scale supersymmetric case [96]. In fact, even the 126 and 10 of Higgs scenario would require complex 10 just 
to meet the charged fermion mass constraints [83] in the non-supersymmetric case. This finally leaves, as the only 
viable possibility, the most general scenario with the 10-, 120- and 126-dimensional representations as the one that 
could accommodate constraints generated by the A phenomenology in the 50(10) framework. The relevant relation, 
in that scenario, reads 

e],DlM'^^^ - M^'^^ElD*j, - -44ri26i?^«i26 + 4E],Y,2oD*j,v';2„. (83) 

Clearly, E]^Yi2gD'^ {E]^Yi2oD'^) would be proportional to Y for A originating from 126 (120). In both cases there 
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are more than enough parameters to accommodate required form of Y. Note, however, that our conservative estimate 
for the upper bound on V45 as shown in Fig. 7 should still be applicable on either z)i26 or w"2o- For example, if we 
identify E'j^Yi2gD^ (£']jYi2o-DJj) with Y it is clear that the left-hand side of Eq. (83) cannot be dominated by the term 
proportional to ^"20 (wi26)- If the opposite was true, we would obtain £^]jFi26-C'Jj ~ i?]jYi2o£'|j which certainly cannot 
hold as Yi2e is symmetric and Y120 is antisymmetric. To conclude, the only viable candidate that can accommodate 
Y is the 50(10) scenario with the 10, 120 and 126. 

VI. CONCLUSIONS 

We have investigated the role of a colored weak singlet scalar possibly addressing the tt FBA puzzle in flavor 
changing processes and precision observables of down-quarks and charged leptons. The magnitude of the predicted 
effects is governed by the mass of the scalar (which we normalize to 400 GeV as preferred by the tt phenomenology) , 
and (generic) complex matrix Y acting in quark and lepton flavor-space. Y is the central object of this analysis. 

Virtual contributions of the considered scalar affect many observables and in order to obtain insight into the Y 
structure we have analyzed a plethora of rare quark and lepton processes, some of them well measured, others bounded 
from above. In particular we have considered FCNC and CP violating observables in K and Bd,s meson systems, 
(lepton flavor violating) dileptonic decays of neutral mesons, /i— e conversion in nuclei, anomalous magnetic moments 
of charged leptons, and lepton flavor violating decays of the muon and r lepton. For completeness, we have also 
considered effects in the Z — > 66 decay width. We have properly accounted for SM contributions to the relevant 
observables where needed. 

Then we have performed a global x^ ^t of the Y matrix elements and found an excellent agreement with all the 
considered constraints. We have confirmed that the couplings to electrons are strongly suppressed. The most salient 
finding is the explanation of the anomalous magnetic moment of muon, which requires the muon coupling to a single 
generation down-quark to be of order one. Combined with LFV B and t decay constraints, this leads to strong limits 
on the tau lepton couplings to down quarks which in turn exclude the possibility [52] to simultaneously explain the 
measured large CP-violating mixing phase in the Bs sector or a large enhancement of absorptive mixing amplitude 
ri2s in this model. Even in absence of the a^ constraint, the B^ — B^ mass difference measurements alone constrain 
the relevant leptoquark couplings and a large new absorptive contribution in Bg — Bg mixing cannot be generated. 
Using a value of \Vub\ preferred by the measured branching fraction of i? — )■ tj/ we find that this model can modify 
the Bd ~ Bd mixing amplitude sufficiently to remove the tension between the two observables. However, in this case 
the anomalous magnetic moment of the muon cannot be explained. 

We have systematically implemented all the phenomenological constraints in a class of SU{5) models where all 
the fermion masses are generated at the tree-level to find out that the explanation of the a^ anomaly requires the 
vacuum expectation value of the 45-dimensional representation to be of the order of 10""'^ GeV. This result implies 
that the up-quark couplings, in this setup, are symmetric in nature. This in turn makes predictions for certain partial 
proton decay lifetimes very accurate. We have also shown that the symmetric scenario for the Yukawa couplings in 
the down-quark and charged lepton case is not compatible with the constraints due to the presence of light A and 
discussed implications for the S'O(IO) type of unification. The simplest of possible realizations of both 5*0(10) and 
SU{5) with the symmetric Yukawa sector, that could accommodate observed fermion masses, are shown not to be 
viable unless A is heavy enough not to play any role in low-energy phenomenology. 

We conclude by noting, that the couplings of the leptoquark in question to the matter fields, in the physical basis, 
are always dominated by just one of the entries of the second row of matrix Y. That entry is at least two orders 
of magnitude larger than any other entry. This property puts this particular leptoquark effectively in the so-called 
second generation category. Moreover, as it does not couple to neutrinos, the bound extracted from the recent LHC 
data for the second-generation leptoquarks [95] is truly applicable in this case and reads mA ^ 380 GeV, accounting 
for the reduced A -^ fij branching ratio of order B > 0.7 due to the presence of the A — )■ tj decay channel [2]. This 
and the upper bound on its mass — m^ < 560 GeV — that originates from simple perturbativity arguments thus 
place it in a very narrow window of discovery. 
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Appendix A: One loop contributions of A to Ri, 



We are working in the massless limit m^ = nib = and in d = 4 + e dimensions to regularize UV divergence. The 
first two diagrams in Fig. 3 give the following contribution to the 1-particle irreducible (IPI) amphtude 



9x| - 2xz + log(xz) (3x1 -6xz + 6 log(l + xz)) (Al) 



3e 6a;| 



(47rj^ 
-8/i + 4/2(2xz - 4) + 6Li2(-xz)) - i7r(34 -6xz + 6 log(l + xz)) 



Ai 
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where C^ = r7i^/(47r)^/^r(l — e/2) and /i, /2 are auxiliary functions defined as 



/i = 4 arctan 



+2Re 



/2 = 2W arctan 



[u^(^^ + yxz{^-xz))^U2{^{i'Xz)-\{l^xzWxz{^~xz))}. 



(\l^-lz) 



(A2) 



(A3) 



In addition to graphs in Fig. 3 there are one loop contributions of A to 6-quark self-energy, corresponding to on-shell 
field renormalization of the 6-quark field 






2 1 

^ + 2 



(A4) 



Combining the tree-level SM with IPf diagrams of A and the field strength renormalization we obtain the UV-finite 
amplitude 



A = Zh{A"''° + A^''''') = igz [{g^R + SgB)AB. + glAL 
where the change of right-handed coupling, Sgu, is given in Eq. (6f ). 



(A5) 



